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We explore the physics of a gas of particles interacting with a condensate that spontaneously 
breaks Lorentz invariance. The equation of state of this gas varies from 1/3 to less than —1 and 
can lead to the observed cosmic acceleration without requiring a vacuum energy. The particles are 
always stable. In our particular class of models these particles are fermions with a chiral coupling 
to the condensate. They may behave as relativistic matter at early times, produce a brief period 
where they dominate the expansion with w < today, and behave as matter at late time. There are 
no small parameters in our models, which generically lead to dark energy clustering and, depending 
on the choice of parameters, smoothing of small scale power. 



I. INTRODUCTION 

A wide range of recent observations (see, e.g., Refs. [| 
and 0) suggest a "dark energy" component that has an 
equation of state close to —1 and dominates the recent 
expansion history of the universe. The existence of the 
dark energy is a profound "fine tuning" problem for par- 
ticle physics 3] . Most solutions that can be expressed in 
terms of a Lagrangian invoke the presence of a non-zero 
vacuum energy 0,[5| ■ 

Yet many physicists in the string theory community 
doubt it is possible to have any vacuum energy at all -6| . 
Braneworld theories that invoke the presence of extra 
dimensions have other difficulties in producing a non-zero 
but small cosmological constant |7| . 

We present an alternative mechanism for generating 
an accelerated expansion and negative equation of state 
that does not rely on a vacuum energy, and where no fine 
tuning is required. In our model, a fundamental field 
spontaneously breaks Lorentz invariance by developing 
a non-zero vector expectation value that thus "picks" a 
particular direction in spacetime. Just as in the case 
of the Higgs field in the electroweak theory the vacuum 
energy of the condensate is set to zero. 

The negative equation of state is produced by a gas of 
fermions that interacts with this Lorentz-violating con- 
densate. At high energies, the fermion gas behaves like 
ordinary radiation; however, at temperatures compara- 
ble to the mass of the fermion, this gas behaves as if it 
has a negative equation of state. 

We present the physics of particle dark energy in two 
parts. In Sec. [HI we derive the equation of motion for 
a fermion, ip, coupled to a constant, time-like vector 6 M , 
in curved spacetime. We then show how such a situation 
may arise in a fundamentally Lorentz invariant theory by 
the process of spontaneous symmetry breaking. We then 
show how a gas of such coupled ip particles may naturally 
lead to a negative equation of state, and possibly w < — 1. 

In the remainder of the paper, Sees. II I II to IVII we 



discuss realistic cosmological scenarios in which particle 
dark energy can explain the late-time acceleration. 



II. LORENTZ-VIOLATING COUPLINGS 

We first consider the dynamics of a fermion field, ip, 
coupled to a constant vector with dimensions of mass, 
a case that has been considered in flat spacetime by 
Ref. [8| . We take the Lagrangian to be 



£ = ip {i^Dft — m)ip- (ipja^ip) V 



(1) 



Here ip is a fermion of mass m (assumed to be small), 
and the operator represents the covariant derivative 
acting on the spinor. 

We take the standard FRW metric to be 



ds 2 



dt 2 - a 2 (t)dx 2 



(2) 



where a(t) is the scale factor. We will take the vector 
6 M to be time-like and, moreover, to be normal to the 
homogenous space-like hypersurfaces of the FRW time- 
slicing. Thus, we can write 6 M = (6, 0, 0, 0), and the equa- 
tion of motion for ip is 
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where a is da/dt and the term proportional to the Hubble 
parameter, a/a, comes from the spin connection term 
given by the metric in Eq. [21 (see Ref. [9(.) When k 3> H 
we may discard terms proportional to H , thus recovering 
the flat-space equation of motion. 

We write the wavefunction ip as a spinor multiplied by 
an exponential, being careful to allow for a non-standard 
dispersion relation between the time and spatial compo- 
nents of the momentum. Our ansatz is 

ip a (x,t) = u a {u,k)e- l[oJ ^ R)t - aW ^, (4) 
ip a {x,t) = v a {w,h)S u *® t - a ® rt '%, (5) 

where we have chosen to write k as the physical (not 
comoving) momentum, we have left the dependence of u 
on b implicit, and our ansatz is that a(t)k is constant. 
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After manipulation and again discarding terms propor- 
tional to H, one finds (taking the chiral basis, and up to 
a normalization), 



our new fermion coupling, is 



u + 



-k + b + uj + (k) 

m 




v+ = 




(kTbf 



(8) 



where we have chosen the direction of propagation to 
be along the z-axis and the spinor subscript ± refers to 
states of positive and negative helicity. 

The existence of these solutions implies that the physi- 
cal three-momentum of a i/j particle redshifts in the usual 
fashion. The particles do not "decouple" from the rest of 
the expansion and k oc a -1 . This is to be expected from 
Noether's theorem, since spatial translations, and thus 
momentum, are preserved. As we shall see, conservation 
of stress-energy also requires this behavior. 

Eq. ^ describes a manifestly non-covariant theory be- 
cause 6 M is fixed for all time. However, it is possible to ar- 
rive at such a Lagrangian as a low-energy approximation 
to a theory that develops a non-zero, vector-valued vac- 
uum expectation value. This is a "spontaneous" breaking 
of local Lorentz invariance, analogous to the way in which 
the Higgs field spontaneously breaks the internal symme- 
try SU(3). Early in the universe, at high temperatures, 
we expect both symmetries, Lorentz and SU(3), to be 
restored. 

To break Lorentz invariance spontaneously, one might 
imagine, for example, that b^ is a vector field with a 
standard kinetic term, the above coupling to the field tf>, 
and a (gauge- violating) Mexican hat potential V(b^b^). 

For the chiral coupling we have chosen in Eq. ^ such 
a model would require, among other things, large di- 
mensionless parameters. We will thus use instead the 
"ghost matter" model |ToL Hl| as a concrete example 
of a Lorentz-violating condensate that can be written 
as a Lorentz-invariant fundamental theory. There are 
many recent models to which this analysis may be ap- 
plied with potentially interesting results {e.g., Ref. |12| . 
which has spacelike condensates), and there is no funda- 
mental reason why our model could not rely on vector- 
or tensor- valued condensates. Importantly, two recent 
works 0, Q have discussed the possibility of sponta- 
neous Lorentz violation via a scalar gradient, and the 
consequent breakdown of Lorentz invariance in the dis- 
persion relation. 

The Lagrangian for the ghost condensate, including 



C = P (X) + iip (f 7 /*D /1 
-^757^)^0+... 



m) if) 



(9) 



Here the field <f> represents the scalar field condensate, 
and d^4>d^(f) is written as X for brevity. The expectation 
of d^4> will be our b^. The term F has dimensions of 
mass and is assumed to be of order the Planck scale. The 
ellipsis represents higher order terms in the field <j) which 
stabilize fluctuations in the <p field; their exact form will 
be unimportant for our investigation. 

The function P is the non-standard kinetic term of 
the ghost matter model, and is chosen to have a mini- 
mum when the scalar field derivative, d^ifi, is non-zero 
and time-like. 

Following Ref. [l^, and by analogy with the Higgs 
mechanism, we take P(X) to have a minimum at some 
non-zero, positive X. We expect this to be of the same 
order as the energy scale of the theory, M, and for sim- 
plicity we take d^cf) at the minimum to be (M 2 , 0,0,0) 
in some frame. The most natural choice is to take this 
frame to be coincident with the CMB. Our parameter b 
is then M 2 /F. 

As we shall see, a negative pressure will arise in this 
system even in the absence of a vacuum energy. In order 
to demonstrate this explicitly, we will take the value of 
P(X) at the minimum to be zero. 

The presence of ^-particles will perturb the field </>. 
The classical equation of motion is 



D M (2FP'(X)d»<j> - (^757^» = 0. 



(10) 



We will consider situations in which there is an equal 
number density of ip particles with positive and with neg- 
ative spin along any axis. In this case, the spatial com- 
ponents of the second term in Eq. 1101 vanish, and no 
spatially varying excitations of the <f> field are produced. 
(The creation of cf> fluctuations by a collection of particles 
with net spin has been considered in the non-relativistic 
limit in Ref. \n&.) 

A simple solution to Eq.^Jcan then be found: that for 
which the argument of vanishes, so that 2FP' {X)d^4> 
is equal to (V'757 /i ?/')- Conceptually, the <f) field relaxes 
to its minimum as the ip particles dilute. 

We can now calculate the stress-energy tensor of the 
full system by the Palatini procedure (see Sec. 8 of 



lull system by tne raiatim pro 
Ref. [ill, and Ref. H3-) We find 



Tl xu = -P(X)g^ + | Ulry^D^i) - [D^i^YH) ■ 

(11) 

Considering the second term, we see that the unusual 
dynamics of the (ft-tp coupling cancel each other, and the 
particle physics energy - the energy conserved in interac- 
tions and that appears in the propagator - and gravita- 
tional energy, T 00 , of the -0 particles coincide. Since this 
is seen only when one considers the complete system, in- 
cluding the contribution from the perturbed <f> field, one 
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FIG. 1: The equation of state as a function of spatial mo- 
mentum for a positive (solid line) and negative (dashed line) 
helicity gas with ^-function distribution, and for a blackbody 
centered on k for a positive helicity gas (dotted line), for the 
case m = 6/10. The blackbody solution is a "smoothed" ver- 
sion of the 5-function case. 



might say that the system conspires to enforce the equiv- 
alence principle. 

If M is at least the electroweak scale, F is the Planck 
scale and m is of order M 2 /F, then today the presence of 
ijj particles changes the value of dgcf) by at most one part 
in 10 60 , and the first term P(X) in Eg . II II contributes an 
energy density 10 60 times less than that of the ip particles 
themselves. This is not a "fine tuning" : this extreme 
factor arises from the fact that the energy density of the 
universe today is much smaller than M 4 . We may thus 
neglect this term in our later analysis. 

We can now compute the equation of state for the tp 
particles as a function of their spatial momentum. We 
consider first a universe where all ^-particles have the 
same momentum k, and are either all positive, or all 
negative, helicity. We find the equation of state for the 
two configurations to be 



w±(k) 



k(kTb) 



3(m 2 + [kTb] 2 )' 



(12) 



which we plot in Fig. \I\for the case m = 6/10. We see 
that for positive helicity particles, the equation of state of 
the universe may become negative; indeed, w may even 
become less than negative one, something that can not 
be achieved with ordinary scalar field condensates. 

It is important to note that the negative equation of 
state does not come from a vacuum energy; the potential 
P(X), as noted above, contributes a negligible amount to 
the energy density of the universe. Instead, the negative 
pressure arises in the gas of particles itself. One way 
to see this is as the conservation of spatial momentum 
requiring that the tp particles redshift, combined with 
the fact that for k <b, the energy per particle increases 
with decreasing k. The requirement that the ip particles 
redshift can be found by requiring that T^ v [Eq. be 
conserved. 



v a,± 



(III; 
dk 



kTb 



^Jm 2 + (k =F b) 2 



(13) 



when k is less than 
k, point in opposite 



and for positive helicity particles, 
6, v g and the spatial momentum, 
directions. Thus another, intuitive way to understand the 
emergence of negative pressure is to note that a particle 
leaving a volume in one direction carries carries with it 
a momentum pointed in the opposite direction. 

A somewhat analogous situation to the one in this pa- 
per can arise in MHD systems, where a magnetosonic 
wave can develop a negative pressure [lj| • Note that the 
behavior of the group velocity in our model is qualita- 
tively different from what happens near an absorbtion 
resonance [Tsj . 

Our discovery of a negative w is related to a previ- 
ous result that found a changing, though not a negative, 
w from an analysis of the dispersion relation for a gas 
of particles in a cosmological scenario involving Lorentz 
violation |2fJ. 

Our model, relying as it does on a particular mech- 
anism for spontaneous Lorentz violation, is relevant, 
though not identical, to a wide class of mechanisms de- 
veloped by numerous other authors. For example, [2l| 
examined a Lorentz-violating system where the vector, 
it M , was a fundamental field of fixed norm and not, as in 
our case, the spatial gradient of a scalar field. 

In that situation, important cosmological bounds ex- 
ist from the coupling of u M to gravity; in particular, the 
possibility of an effectively changing Newton's constant G 
arises; see, e.g., |22| and references therein. Our model's 
(minimal) coupling to gravity, which we have derived 
above, is more complicated. 

Later papers that consider a fixed- norm , including 
[2^| . have further developed the theory, considering the 
nature of perturbations in the u field. In our model, as 
discussed earlier in the section, the high mass scale of the 
<j> field means that late-time production of perturbations 
will be negligible; these perturbations certaintly exist, 
however, and their nature is left for future investigation. 

In Fig. |21 we plot v g as a function of k. As the system 
is non-dissipative, the sound speed, c s , is just v g /V3. 
Comparing this with Fig.^ we see that c s — > as w — > 0; 
when the ip particles have the same equation of state as 
ordinary dark matter they also have the same clustering 
properties in the linear regime. This suggests that it may 
be possible for the ip particles to serve some of the roles 
of both dark energy and dark matter, a possibility we 
examine below in Sec. IIVI 

Generically, we expect the ip particles to be in thermal 
equilibrium at some very early time. At some later time, 
the ip particles decouple, and their distribution then red- 
shifts. The decoupling could happen either early, i.e., 
when the ip particles have T 3> m, or late, when T « m. 
We consider the two cases separately, now in the context 
of a realistic cosmology. 
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FIG. 2: The group velocity, v g , for positive helicity (solid 
line) and negative helicity (dotted line) particles, for the case 
m = b/W. 



III. EARLY DECOUPLING 

Let us first consider the case where the ip particles 
decouple while still relativistic. In this scenario, a Fermi- 
Dirac distribution simply redshifts without changing its 
shape. 

From an examination of Eq. 1121 we see that for the 
equation of state to become negative, b must larger than 
m. This restricts the number of e-folds in k over which 
the particles may be distributed; as b becomes larger, w 
can be more negative, but over a smaller Ak/k. 

The width in k over which w is significantly negative 
is fcfmish/fcstart ~ (b + m)/(b— m) < 1.2. Since the Fermi- 
Dirac distribution at large temperatures is significantly 
wider than this, the possibility of a very negative w, and 
universal acceleration, is hidden. This is shown in Fig.^ 
where the equation of state as a function of scale fac- 
tor for a redshifting Fermi-Dirac distribution of positive 
helicity particles is plotted as the dotted line. 

We find a brief period of mild dark energy behavior be- 
tween the radiation and matter-like eras. Including both 
positive and negative helicity states reduces the positive 
and negative excursions of the equation of state further. 

It may be checked that if the ip particles are subdom- 
inant to cold dark matter at small scale factors, when 
they have w rs 1/3, they are always subdominant. Thus 
there is no way to produce a significant period of mat- 
ter domination before the onset of a negative equation of 
state. Early decoupling, then, is not a promising model 
for the late-time acceleration, although it may have some 
applicability in the early universe. 



IV. LATE DECOUPLING 

We now consider the case opposite to early decoupling. 
We assume that, due to a significant self-interaction term, 
the ip particles decouple only after they become non- 
relativistic. We approximate this situation by assuming 
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FIG. 3: A series of "snapshots" of the particle distribution, 
n(k) (solid lines) for the case m = 1, b = 10. The dashed 
line is the dispersion relation, ui(k). Reading the solid lines 
chronologically from right to left, this figure shows how a 
high-temperature Fermi-Dirac distribution cools, and parti- 
cles clump in the well of the dispersion relation at k = b and 
then redshift after decoupling. 

that before the interval between interactions is equal to 
l/H, the particles remain in exact thermal equilibrium, 
at a temperature T corresponding to the energy density 
of particles. When the interaction rate is equal to l/H, 
we assume that the particles gradually drop out of equi- 
librium over a finite range of scale factor. 

We find that before decoupling the positive helicity 
particles form a narrowing clump in the "well" of the 
dispersion relation at k — b where w = 0. The population 
of negative helicity particles is strongly suppressed. 

While the ip particles remain in thermal equilibrium 
the universe appears to be undergoing a standard tran- 
sition from radiation to matter domination. When the 
particles drop out of equilibrium, the very narrow dis- 
tribution at k — b redshifts, and we recover a strongly 
negative equation of state for roughly one e-fold. At very 
large a, the universe is matter-dominated. 

Fig. illustrates this with three overlaid "snapshots" 
of the ip particle distribution. The figure is to be read 
right to left in chronological order. The first snapshot 
is well before decoupling when T w 100m and the par- 
ticles behave approximately like radiation. This is the 
same as the early decoupling case, since relativistic par- 
ticles maintain a thermal profile even after decoupling. 
This distribution is nearly indistinguishable from that of 
photons in thermal equilibrium at the same temperature. 

The second snapshot is before decoupling but after 
the universe has cooled further; here T w 0.25m and 
w ~ 0. The evolution has diverged from the early decou- 
pling case, and the particle distribution forms a narrow- 
ing clump at the bottom of the well. 

The third snapshot is after the system has decoupled. 
The distribution has redshifted towards k = 0, and now 
w < 0. We take decoupling to happen over a finite range 
of scale factors: different particles leave the bottom of 
the well at slightly different times, so the distribution is 
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slightly wider than just before decoupling. 

We can thus have an arbitrarily long period of matter 
domination between the radiation and dark energy dom- 
inated eras. As before, the dark energy era lasts only 
a brief period of time, and now appears soon after the 
interaction time becomes of order l/H. In contrast to 
the early decoupling case above, however, a previously 
subdominant tp-Emd component can grow to become the 
dominant contribution to the stress-energy. 

Let us now consider the plausibility of a late decoupling 
around z « 0, to explain the contemporary acceleration. 
We require that 

noav ~ Hq (14) 

where Uq is the number density of particles, of order 
Pa/m, a is the cross section, which we take equal to 1/m 2 , 
and v is the average particle speed. 

It remains to determine v. To a first approximation, 
we can take v approximately equal to c at all times. This 
seems counterintuitive, because if a particle is kicked into 
a state k « b, its group velocity is very small. The par- 
ticle may drop out of equilibrium and start to redshift. 

However, the unusual nature of the group velocity im- 
plies that a redshifting particle at k = b finds its veloc- 
ity increasing. The steepness of the group velocity curve 
around this point means that a particle need only redshift 
by a factor of 1 + 0.1m /b in order for its velocity to rise 
from zero to « 0.1c. After redshifting this small amount, 
it can start to reestablish the equilibrium suggested by 
setting v to unity in Eg. ITU above. 

Setting v above to unity, then, we can find the maxi- 
mum mass of the tp particle such that decoupling happens 
around to today. Nucleosynthesis gives us the constraint 
that, if the tp particles are to make up a non-negligible 
portion of the energy density today, they must enter the 
w = phase before T ss 1 MeV. We find, to an order of 
magnitude, m < 100 MeV if the particles are to decou- 
ple today, which easily satisfies this bound. Setting the 
fraction of tp particles to fi^, ^ 1, or changing the cross- 
section by a factor a 2 lowers this mass estimate only by 
a factor a 2 / 3 ^ 3 . 

We thus consider a "^CDM" model, a flat universe 
with similar energy densities of tp matter and some cold 
dark matter particle. Approximating the tp decoupling 
as happening over some finite range in scale factor adeo 
we plot in Fig. 0] the equation of state of the universe for 
$CDM. with ACDM plotted for comparison. Following 
the discussion above, we take decoupling to happen over 
a range of redshifts; the particular choice here is 0.20 < 
z < 0.62, and we take fi^, before decoupling to be 0.08. 

We have chosen parameters so that both models have 
the same Hq, d\ ss and deceleration parameter qo; there 
are many different combination of decoupling redshift, 
O^, and b that can satisfy this requirement. Making the 
decoupling process take longer widens the "well" of the 
effective equation of state in Fig. Q] 

The simple analysis above suggests that it is possible 
to describe the late-time acceleration of the universe with 
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FIG. 4: The equation of state of the entire universe (i.e., 
including all components) as a function of scale factor in a 
'l/CDM scenario (solid line.) Here we have taken the decou- 
pling to happen between 0.20 < z < 0.62, we have taken the 
Qy,^ = 0.08 at z > 1, and b = M 2 /F = 400m. The d LSS , 
h and qo for this model agree with that of a fiducial ACDM 
model, plotted for reference as a dashed line. 



a very simple model that relies on a late-time decoupling 
of the tp particles. There is no fine-tuning required: the 
tp particles must be subdominant to satisfy large scale 
structure constraints, but their contribution to the den- 
sity before decoupling may be between 1% and 10% of 
the CDM. If our tp particles are to have mass near 100 
MeV, then if F is the Planck scale, the intermediate scale 
M is ~ 10 9 GeV. 

The fraction of energy density contained in tp parti- 
cles at high redshift is non-negligible, and depending on 
choice of parameters can be of order tens of percent. At 
high redshifts, the ip particles act like ordinary dark mat- 
ter, and can cluster. Closer to today, during the contem- 
porary acceleration, they have a near-unity group ve- 
locity. In many ways, then, the phenomenology of our 
model is similar to a decaying dark matter scenario |2-l| , 
where some fraction of the dark matter turns into ultra- 
relativistic particles, smoothing out small scale power. 

It is this freestreaming that makes it difficult for the 
tp particles to be both dark energy and all of the dark 
matter, although more sophisticated models, with a more 
detailed treatment of the unusual decoupling process and 
the evolution of perturbations, may allow us to unify 
these two mysterious components. We consider this an 
open question; we discuss the clustering properties of this 
model in greater detail in Sec. IVII below. 



V. DECAYING DARK MATTER 

We believe that the late-time decoupling scenario is 
the most promising particle dark energy model of the 
kinds discussed in this paper. However, it is instructive 
to consider other scenarios, in part so as to understand 
the particle physics aspects of the tp fluid. 



6 







\\ 

v\ 
\ 

\ 

\ 


\ 

\ 

\ 



-2-10 1 2 3 
log 10 a 

FIG. 5: The equation of state of the entire universe as a 
function of scale factor in a D^CDM scenario (solid line.) 
Here we have taken r = 2io and b — 45m. The dLSS and h 
for both models agree with that of the ACDM model, plotted 
for reference as a dashed line. 



Our second toy model is for the tp particles to be pro- 
duced during the late-time decay of a separate dark mat- 
ter particle. Our so-called "D'FCDM" model is a modi- 
fied version of the decaying CDM models first proposed 
in Ref. |24|, where in D^PCDM the decay products are tp 
particles. 

Let us assume that the dark matter particles are de- 
scribed by a single scalar field \ °f mass m x that may 
decay via a -0-number conserving interaction with an ex- 
ponential decay time r we take to be of order 1/Hq. We 
assume the x particles are cold. 

If the mass of the x particle, m x , is greater than 
V m 2 + b 2 , the x particle decay products will include both 
positive and negative helicity tp states. However, if m x 
is less than s/rr? + b 2 , decay into negative helicity states 
is energetically forbidden. 

Let us consider this case and take, for simplicity, m x 
equal to m, the mass of the tp particle. Then, the % 
decay products end up in a narrow distribution of positive 
helicity tp particle states at k = b. As the tp particles are 
produced, they redshift to k < b and, as above, develop 
negative pressure. 

We approximate the universe as beginning with only x 
particles; the tp population is built up as the x particles 
decay with e-folding time r. Fig. [5] shows the equation 
of state of the universe, as a function of redshift, for a 
model where we have taken r to be of order the age of the 
universe, to, and have chosen the mass scale, b, such that 
the distance to last scattering, and the Hubble constant 
today, are the same as a fiducial ACDM, which we plot 
for comparison. 

Though both models, D1-CDM and ACDM, in Fig. 
have the same c?lss and h, they give different values 
for the deceleration parameter, q . For D'J'CDM, q w 
-0.05, while for ACDM q sa -0.55. 

We find thus that this decaying dark matter model 
naturally produces a period of accelerated expansion in 



the low-redshift universe even if only a small fraction of 
dark matter has decayed. This is because the equation of 
state for the tp particle population is very negative just 
around k ss b; the tp particles grow rapidly in comoving 
energy as soon as they are produced. 

The position of the turnover depends both on the rel- 
ative values of b and m, and on the dark matter lifetime. 
One can keep the position of the turnover fixed by in- 
creasing both b/m and the dark matter lifetime, which 
the overall shape remains unchanged. The equation of 
state returns to zero after most of the dark matter has 
decayed and the majority of tp particles have redshifted 
to k < b. 

There are three mass scales in our theory: M, F and 
to, where to is supposed to be near M 2 /F. If we take 
F to be the Planck scale, then if M is the GUT scale, 
the dark matter mass is greater than 10 12 GeV. Alter- 
natively, if our dark matter is an Axion condensate with 
to ~ 10~ 5 eV, M may be as low as 100 TeV, roughly the 
electroweak scale. 



VI. EVOLUTION OF PERTURBATIONS 

The evolution of perturbations in particle dark energy 
models remains to be studied; however, broad statements 
may be made. Both *CDM and D*CDM make definite 
predicitions about the growth of large scale structure. 

First, both predict the signatures of dark energy clus- 
tering. In the 'FCDM model, the sound speed, related 
to v g , the group velocity, goes to zero with w just be- 
fore the acceleration, so that the tp fluid perturbations 
may grow like ordinary matter before their equation of 
state becomes negative. In D^CDM, the progenitors of 
the tp particles may cluster like ordinary dark matter. 
The reader is cautioned that the behavior of clustering 
dark energ y in the non-linear regime requires a separate 
treatment |25|. 

Whether and how dark energy clusters remains an ur- 
gent question both for observation and theory pol Wf\ ; 
tantalizingly, we may have seen evidence for it in the 
deficit of law-t power in the CMB p^ . 

Second, both models may have direct relevance to the 
possible excess small scale structure problem (see, e.g., 
Ref. (2(|.) In 'J'CDM, the tp particles may make up a 
substantial fraction (of order tens of percent) of the en- 
ergy density of the universe. When they decouple, their 
group velocity rises rapidly to unity, and they may stream 
out of overdensities on small scales, smoothing out small 
scale power. 

This is just as in the first decaying dark matter mod- 
els [24|. D'FCDM, of course, actually is a model with 
decay: there dark matter turns into ultra-relativistic par- 
ticles that then develop a negative equation of state. 
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VII. CONCLUSIONS 

We have discussed the phenomenon of Lorentz- 
violating couplings in the context of an FRW curved 
spacetime, and have demonstrated how a fermion cou- 
pled to a spontaneous Lorentz-violating condensate can 
produce a strongly negative equation of state. We have 
also presented two models, one requiring a late-time de- 
coupling ("^CDM"), and the other requiring a late-time 
decay ("D^CDM"), that, at least in outline form, can ex- 
plain the contemporary acceleration. In particular, the 
"fCDM model requires no fine-tuning, and involves nat- 
ural cross-sections. 

The rich set of behaviors, and the number of scenar- 
ios that can be constructed even from such a toy model 
as Eq. ED suggests that the question of using different 
Lorentz breaking mechanisms, such as a non-zero expec- 
tation value for a vector or tensor field, deserves further 
investigation. There are also a wide range of possible 
couplings between such a Lorentz-violating expectation 
value and various kinds of matter, including not only 
fermion couplings, but also derivative couplings to scalar 
field particles, and direct couplings to particles of higher 
spin. 

A catalogue of Lorentz violating couplings to the stan- 
dard model has been investigated in Ref. 30]; here, 
since we are coupling to the dark sector, we may 
use much larger couplings without violating Lorentz- 
violation bounds from atomic or accelerator physics. A 
first indicator of whether or not a system will generate 
a negative equation of state in any model of interest will 
be a turn-over in the particle's group velocity, Eq. 1131 

Lorentz violation has long been a subject of study 
in the particle and atomic physics community; we are 
only just beginning to discover its cosmological implica- 



tions |3l|. Many recent astrophysical studies have fo- 
cused on high energy signatures, such as the GZK cut- 
off [32| , and in the context of loop quantum gravity or 
non-commutative geometry |33L l34l l35| . Here, in con- 
trast, we have demonstrated how Lorentz-violating cou- 
plings can lead to negative pressures in the context of a 
spontaneous breaking of Lorentz invariance, and we have 
shown how such violation may have astrophysical impli- 
cations at low energies. 

Spontaneous Lorentz violation has long been consid- 
ered as a feature of string theory for the bosonic string 
and the superstring (3(1 H3, Hg . It may also be a generic 
aspect of braneworld models, where Lorentz violation oc- 
curs in the bulk and may be seen on the brane as well |39j| . 
Finally, the role of the universal time-like vector could 
even be played by the same mechanism underlying the 
thermodynamic "arrow of time." 

Because of the generic nature of our discussion, we 
expect the essence of our conclusions to hold true for 
these more elaborate scenarios when a time-like Lorentz- 
violating condensate couples to the dark sector. As this 
paper has shown, there is a rich set of phenomena to be 
discovered that may potentially resolve long-standing 
theoretical and observational questions. 
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